


pd<-paste(wd, "/Figures/Figures_Output/",sep="")


china<-read.csv("china_original.csv")


chinacf<-read.csv("china_cfactual065.csv")

chinacf4<-read.csv("china_cfactual11.csv")




pdf(file=paste(pd,"chinese_beliefs.pdf",sep=""), 6, 6)




par(mar=c(2,4.25,2,1))

plot(1980:2010,china[30:60],lty=1,lwd=2,ylim=c(-.035,.066),axes=F,xlim=c(1980,2010),type="n",xlab="",ylab="")
mtext(expression(beta[CHN]^{D==1}-beta[CHN]^{D==0}),2,line=2.5,cex=1.25)

grid()

lines(1980:2010,china[30:60],lty=1,lwd=2)


lines(2001:2010,chinacf,lwd=2,lty=2)


lines(2001:2010,chinacf4,lwd=2,lty=3)

axis(1);axis(2,seq(-.04,.06,by=.02),labels=seq(-4,6,by=2))


mtext("Chinese Beliefs",3,line=0,font=2,cex=1.5)

text(2001,.05275,"11.0%",cex=.75)

text(2005.25,.049,"6.5%",cex=.75)



dev.off()







###################################

dprk<-read.csv("dprk_original.csv")[,1:60]


dprkcf<-read.csv("dprk_cfactual065.csv")[,1:10]

dprkcf4<-read.csv("dprk_cfactual11.csv")[,1:10]



pdf(file=paste(pd,"dprk_beliefs.pdf",sep=""), 6, 6)




par(mar=c(2,4.25,2,1))

plot(1980:2010,dprk[30:60],lty=1,lwd=2,ylim=c(-.035,.066),axes=F,xlim=c(1980,2010),type="n",xlab="",ylab="")
mtext(expression(beta[PRK]^{D==1}-beta[PRK]^{D==0}),2,line=2.5,cex=1.25)

grid()

lines(1980:2010,dprk[30:60],lty=1,lwd=2)


lines(2001:2010,dprkcf,lwd=2,lty=2)


lines(2001:2010,dprkcf4,lwd=2,lty=3)


axis(1);axis(2,seq(-.04,.06,by=.02),label=seq(-4,6,by=2))


mtext("North Korean Beliefs",3,line=0,font=2,cex=1.5)

dev.off()
